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ABSTRACT 

A regression technique has been developed to forecast 24-hour char ges of the maxi- 
mum winds for weak (maximum winds ^ 65 kt) and strong (maximum winds > 65 kt) 
tropical cyclones by utilizing satellite measured equivalent blackbody temperatures (Tgg) 
around the storm alone and together with the changes in maximum winds during the pre- 
ceding 24 hours and the current maximum wi .ds. Independent testing of these regression 
equations showed that the mean errors made by the equations are lower than the errors in 
forecasts made by the persistence techniques. 
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PREDICTING TROPICAL CYCLONE INTENSITY USING 
SATELLITE MEASURED EQUIVALENT BLACKBODY 
TEMPERATURES OF CLOUD TOPS 

1.0 INTRODUCnON 

Hurricane caused damages in the United States average over $600 million per year. 
Damages from tropical cyclones arc even much greater in many other countries on the 
>vcstcrn borders of the tropical oceans when expressed ns a percentage of the gross national 
product. Forecasts of hurricane occurrence arc, therefore, of great importance. 

Observations and forecasts of hurricane wind speeds arc also very important because 
damages caused by hurricanes vary exponentially with the maximum wind speeds. WhUc 
the force of the wind varies with the square of the speed, some of the historical surveys 

•V ’• 

of total storm damage suggest that the latter varies with a higher power of the wind speed, 
i.e., 

D = K V“ 

where D is the total damage caused by the storm, K is a constant, V is the maximum wind 
speed and n is some number between 2 and 5 (Howard, ct al., 1972). This relationship 
emphasizes the importance of knowing the intcns&y of a tropical cyclone. 

In spite of the need for knowledge of tropical cyclone intensity by the hurricane fore- 
cast services, aircraft reconnaissance of tropical storms is being reduced to save money. 
Efforts have been increased in recent years, therefore, to use satellite data to observe and 
predict the intensity of tropical cyclones. Results from these efforts have been encouraging 
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iind they keep expanding us satellite data of improved quality become more readily avail- 
able with each new satellite series. 

Satellite measured equivalent blackbody temperatures (Tqd) of cloud tops in tropical 
cyclones should contain useful information about storm intensity and expected changes of 
intensity. Latent heat released when the warm moist tropical air ascends in major cumulus 
towers of hurricanes is the primary fuel for the storms (Dunn and Miller, 1960), and its 
availability is indicated by measurements of the amount and vigor of the convection within 
the cyclone which can be deduced from temperature measurements. 

The hurricane is a prolific producer of clouds. The convective towers build far into 
the troposphere and sometimes penetrate the tropopause, thus producing very cold cloud 
tops. Tlie high level shearing and outward spiraling winds spread the cold cirrus over a 
large area beyond the region of most active convection. This air subsides as it spirals away 
from the storm center causing the cirrus to begin dissipating as the air warms adiabatically. 
These spiralling subsiding effects arc easily observable in satellite imagery and can be quan- 
tified through measurement of cloud top temperatures. Thus, areal distribution of the T 3 B 
provide information showing the extent and strength of the convection wliich serve as in- 
dices of the latent heat released and indicate the extent that the clouds of the storm are 
organized into patterns. 

The latent heat is ultimately converted to the kinetic energy which causes the extreme 
winds of the tropical cyclones (Riehl, 1954). For this to take place, however, there is a 
complex process involving among other things conversion of tlic heat to potential and 
available potential energy. Finally, the kinetic energy has to be concentrated by the flow 
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patterns (mo relatively narrow bands for the storm to become truly destructive. All these 
processes take time and there should be a lag between changes in convective activity and 
changes of maximum winds in the storm. 

The results of theoretical^numcrical model experiments simulating development and 
maintenance of tropical i' /clones support this reasoning and suggest that maximum vertical 
motion, that is, maximum convection, prjjcedes the maximum winds by one to tlircc days 
(e.g,, Rosenthal, 1978; Kurihara and Tulcya, 1974). Riehl (1954) and Rosenthal (1978), 
have also emphasized that the convection needs to be organized by some larger scale sys* 
tern into a suitable pattern (c.g., spiral bands and eycwall before rapid intensification of 
the tropical cyclone takes place). 

Dvorak and earlier investigators at the National Environmental Satellite Service have 
developed techniques to use satellite imagery to identify the present intensity of the trop- 
ical cyclone and to suggest future changes of the intensity (Dvorak, 1973; Hubert, et al,, 
1969). Wliile these teclmiqucs have shown skill and the latest Dvorak teclmique is in wide- 
spread use in the tropics v/orldwide, it still involves considerable subjectivity especially in 
the forecasting of storm intensity. Dvorr.k (1973) utilized the degree of pattern organi- 
zation to identify the current storm intensity from satellite imagery. He found that the 
size of the central dense overcast of cirrus and the degree to which the spiral bands of 
convective clouds encircled the storm center to be important factors. 

ir 

Based on the heuristic reasoning just presented, results from the theoretical experi- 
ments, results of using Dvorak’s technique under operational conditions, and other research, 
the authors developed a hypothesis to be tested by the experiments reported in this paper. 
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It says: (I) The Tyjj of the tropical cyclone cloud tops provide a measure of the convection 
and an index of the latent licat released I Dr eventual conversion into kinetic enerBy; (2) The 
Tpu areal distribution serves as an Index of the organization of the storm’s convective activ- 
ity; and (3) The lower (higher) the mean Tgu of the cloud tops over a moderate sized area, 
the stronger (weaker) and more (less) persistent is tlic convection and the more likely that 
tiie maximum winds in the storm will increase (decrease) with time, 

2.0 THE DATA AND THE ANALYSIS 

The Tjjjj for a number of tropical cyclones were analyzed using data from the Western 
Atlantic for 1969 and the Western Pacific for 1970, 1973 and 1974. TheTgj) in 1969 were 
measured by the MRIR sensor on NIMBUS 3 with a spatial resolution of 55km (at nadir). 
Those in 1970, 1973 and 1974 were measured by the THIR sensor on NIMBUS IV and V 
with a spatial resolution of 8 km at the subpoint. These data were analyzed using the 
scheme illustrated in Figure 1 to get a measure of the intensity, expanse and organization 
of the storm. The concentric circles are 111 km apart and the rings they bound are num- 
bered outward from I to 12. The mean temperature was computed for each ring with the 
center of tiie diagram coinciding with the center of the storm. In addition, the mean tem- 
perature was computed for each octant of each ring (hereafter referred to as a sector) with 
the top of the diagram being oriented both towards the north and also along the direction 
of motion of the storm. To get a further measure of how well the convective towers were 
distributed symmetrically and concentrated about the storm center, the standard deviation 
of the mean sector temperatures were computed for rings 1 through 5 and for various com- 
binations of rings. With these detailed data it is feasible to study the expanse and also the 
organization of the storm as well as tiie intensity of the convection. 
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3.0 Ti,u AND STORM INTENSITY 

Tiic first tests made in this investleation were with 1969 Allmitlc tropical cyclones 
and lielpcd evaluate how well the temperatures demonstrated an index of convection. The 
mean data for the rings composited for 16 hurricanes are compared with similar data for 
19 storms of less than hurricane intcu'iUy in Figure 2, For rings 1-4 (Fig. 1) the mean 
temperatures were 7 to 10®C lower in the hurricanes than in tlic weaker storms, but in 
rings 6-10 tlic luirricanes had higher temperatures. Tliat is, the temperatures imply tliat 
hurricanes have stronger convection near the core and stronger subsidence in the environ- 
ment .surrounding tiie storm. Both the convective and subsidence areas have frequently 
been obsen’cd by aircraft reconnaissance and in satellite Imagery (Slienk and Rodgers, 
1978). The subsidence dissipates many of the clouds at distances greater than 650km 
from the center and tiuis causes the higher mean temperatures. 

A simitar test was made witli the 1970 Western Pacific tropical cyclones. Data were 
composited for 3 groups of slormi those of less than typhoon intensity (15 cases), ty- 
phoons with maximum winds less than 100 knots (13 cases), and typhoons with maximum 
winds equal or greater than 100 knots (14 cases). The data are graplicd in Figure 3, and 
the differences are plotted In, the insert. The comparison between the weak storms and 
those with maximum winds greater than 99 knots is similar to that of the 1969 Atlantic 
storms except for ring 1 where the typhoons were warmer. Tills reflects tlie fact that 
many of tlie typhoons had large cloud free eyes and tlie THIR instrument used in 1970 
on NIMBUS IV had sufficient resolution to measure high temperatures over smaller areas 
than tlie MRIR sensor used in 1969 on NIMBUS HI. If values for ring 1 are ignored, in- 
tense storms are again colder in the inner four rings and warmer at greater radii. (Broken 
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Figure I. Grid used in analyzing the tem(>erature data. Tlie con- 
centric rings are spaced 1 1 1 kni apart. The center of this grid 
coincides with the center of the tropical cyclone. 
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mNGS.-.DISTANCE ROM STORM CENTER 

Figure 3. Comparison of mean temperatures of cloud tops around I <>70 trop- 
ical cyclones of the Western North Pacific: IS storms with maximum winds 
less than 65 knots (broken line); 13 typhoons with maximum winds less than 
100 knots (long-short-long dashes); and 14 typhoons with maximum winds 
equal or greater than 100 knots (solid line). The insert contains graphs of the 
differences in temperatures of the latter and the two weaker categories. The 
rings arc illustrated in Figure I. 
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line labeled 3-1 in the insert.) Tlie comparison between the typhoons of moderate inten- 
sity with the very intense t^yphoons (solid line, 3-2. of insert), however, gave contradictory 
information. An examination of the individual cases revealed that the sample of moderate 
intensity storms were biased toward storms that were intensifying and the more intense 
typhoons were biased toward mature storms thjit wore changing slowly or were weakening. 
This is especially significant because it further suggests that the mean temperatures are an 
index of the rate of change of storm intensity. 

4.0 FORECASTING STORM INTENSITY 

In Figure 4 the 1970 data for storms south of 30'*N arc stratified according to the 
change of intensity during the succeeding 24 hours. The four categories used were inten- 
sifying (maximum winds increasing at least 10 knots - 12 cases), weakening (maximum 
winds decreasing at least 10 knots - 4 cases), little change (includes all storms which 
eventually reached hurricane intensity where the change during the next 24 hours was less 
than 10 knots - 6 cases), and a group of tropical storms which never intensified to the 
typhoon stage before they finaily dissipated (5 cases). We can note that the storms with 
the greater rate of intensification arc associated with the lower mean temperatures in ali 
rings out through 9. The intensifying storms are about 18“G colder in rings 2-4 than the 
storms that never reached typhoon intensity, ibout 10‘^C colder in rings 1 through 8 than 
the weakening storms, and 5“ to 10'*C colder than the storms changing slowly in intensity. 

Figure S illustrates a time-lag between the Tg^ of the cloud tops and the changes of 
the maxil aum winds in two typhoons. The temperature scale is inverted to show the lower' 
temperaCures at the top of the diagram. To interpret the graphs one might consider the 
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TEMPERATURES *K 



Figure 4. Same as Figure 3 except that the data are stratified according to 
change of intensity during next 24 liours: Maximum winds increasing 10 or 
more knots (1), maximum winds changing less than 10 knots (S), maximum 
winds decreasing 10 or more knots (W), and storms which never reached 
hurricane intensify (T.S.). Only storms located south of 30®N were iit- 
cludeJ. The insert shows that the intensifying siorms have much colder 
cloud tops within 8 degrees (888km) of the storm center than the others. 
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Figur;, 5. Temporal changes of mean equivalent blackbody temperatures (nrigs 2, 3 and 4 from 
Figure I) and maximum winds of Typhoons Hope (Sept. 1970) and Billie (August 1970). The 
maximum wind changes lag the temperature changes in both cases. The temperature scale is 
inverted to show the lower temperatures at the top of the graph. 




temperature graph as a crude index of the convective activity. Note that the minimum 
temperature (or maximum of tlje convective index) for Typlioon Billie occurred more than 
2 days earlier than the maximum winds. 

In the case of Typhoon Hope the wind graph maximum lags by more than one day 
(missing Tjjg data make it inypractlcal to determine tlie exact time of the temperature 
minimum). Similar data have been examined for other storms and results suggest that 
changes in tiie maximum winds lag changes in the temperatures by 24 to 36 hours. From 
this we can conclude that tiic Tbb contain predictive information. 

Figure 6 is adapted from a simulation experiment with a theoretical model by-> 
Rosenthal (1978), The vertical velocities ,it 900inb and tlie maximum winds are plotted 
against time. Hero vertictii velocity, rather than temperature as in Figure 5, represents 

convection. Note the similarity in the time iag between maximum convection and the 

i 

maximum winds of the slorm in the two illustrations. |l. 

Analyses of the data and heuristic reasoning suggest that the Tpu will viury with at 
least tlie following: rate of cliange of intensity of the storm, intensity of the tropical 
cyclone, latitude of tlie storm, season of the year, and meafi temperatures tlu-ougli die 
surrounding troposphere. Tlie latter means that relationships may differ for tlie various 
oceans. Because of these and other factors, the data used in preparing the composite 
graphs reproduced in Figures 2, 3 and 4, show considerable scatter. It was necessary 
therefore to consider otiier parameters as well as the Tgjj when developing an objective 
teclmique for forecasting clmnges in tropical cyclone intensity. Even a cursory examin- 
ation of tlie data reveals, for example, that the relationship between Tj)b and future 
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Figure 6. Comparison of Vertical Motion and Central Pressure in a Model Hurricane 
Showing Lag with Time (Rosenthal 1 978) 



storm ‘intensity is different between storms that have been intensifying and storms tiiat 
liavc been weakening in the preceding 12 to 24 hours. It also seems obvious that the 
arc indicators of both the current and future intensities of the storm. These two effects 
need to be separated if a higiUy successful predictive scheme is to be developed. The 
change in maximum wind speeds during the preceding 24 hours and the current maximum 
wind speed helped define these effects.* A screening regression procedure selected the Tpg 
and these latter two parameters in determining the regression equations which gave the 
best results for forecasting maximum winds of tlie tropical cyclones. 


Two regression equations have been developed to predict the changes in maximum 
winds of tropical cyclones during the next 24 hours; the first is for the weaker storms 
where the current maximum winds arc equal or less than 65kts, and the second is for the 
more intense storms. 

V+24w = ‘43.75 - 0.S94T2,3 + 0.389AV.24 + Vq (1) 

j ^24, = 227.86 - 0.76Ti,2.3 + 0.499AY.24 + 0.398Vo (2) 

where 


V+24w (V+ 24 j) are predicted maximum wind speed (knots) for storms whose current 
maximum wind <65 knots (>65 knots). 

:i 

^2,3 (^ 1 ^ 2 , 3 )^*'® between 1 10 and 330km (0 and ;^30km) about the 

storm center [rings 2 and 3 (rings 1, 2, and 3)]. 


*the winds were extracted from the “best track" information published in the Joint Typhoon Warning 
Center Annual Typhoon Report for 1970, 1973, and 1974 (Reference 11, 12, and 13). 
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JiV li the change in inaxinium wind speed of the storm during the preceding 24 
hours (knots). 

it the current maxiinuni wind speed (knots) in the storm. 

These equations were develo|H*d using as dc|X*ndent data 5H cases from the 1 470 
Western Pacific tropical cyclones. Pertinent statistical information is summari/ed in Table 
I and show that errors of forecasts made with the alnive regression equation* are consider- 
ably less than those of techniques using only persistence. These latter techniques have 
been used by forecasters to- many years, and except in s|K*cial situations give about as 


Table I 

Mean Trrors* of T'orecasts (Knots) 


IX'pi'ndcnt Data 
(D>70 Storms) 

Independent Data 
(l‘)73and 1974 Storms) 

lc| nation 

Vo 

n 

12 

t;-nc'2 


n 

i;2 

H-NCJ 

H-P< 

1 

^ 65 

24 

8.6 

-4.0 

-4.2 

33 


-4.2^ 


'> 

> 65 

34 

y.4 

-7.5 

-2.6 

20 

16.1 

-3.1 


1 & 2 

> 25 

58 

«J.l 

-6.0 

-.T2 

53 

1.3.9 

-3.8^ 

D 


I. Mean errors are the means (in knots) of the absolute values of the differences between 
forecasted and observed maximum wind speeds. 


2. I' is the mean error of forecasts made by the indicated regression equation. 

3. N(' is the mec. error of forecasts made by assuming no change in wind speed during the 
forecast periinl. 

4. P is the mean error of forecasts made by assuming the wind speed would change the 
same amount during the next 24 hours as it did during the preceding 24 hours. 

•Differences are statistically significant at the I percent level. 

••Differences are statistically significant at the 5 percent level. 
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good results as any tccimiquc currently used (Gentry, 1973). They are: (1) persistence 
forecast assuming that no change would take place during a forecast period (NC) and (2) 
persistence assuming that the change during the next forecast period would be the sjunc 
as the change during the preceding period (P). 

Tests of Equations (1) and (2) with independent data were made using 56 storm cases 
from years 1973 and 1974. Statistics of the mean errors are summarized in Table I. The 
errors in tlic forecasts made by Equations (1) and (2) are again lower than errors in fore- 
casts made by the persistence techniques. The difference between the errors of Equation 
(1) and those made by forecasting no change, and the difference between errors of fore- 
casts by Equation (2) and forecasts of persistence of change are both significant at the 1 
percent level. The combination of Equations (1) and (2) make forecasts for all the storms 
significantly superior to either of the persistence techniques. 

5.0 FORECASTS USING ONLY SATELLITE INFORMATION 

The equations just evaluated used parameters other than those measured by satellites 
Tlie question may well be asked whether skill can be shown using only the satellite infor- 
mation and without including the change in wind speeds during the preceding 24 hours as 
a predictor. To answer this question, equations were developed to predict the maximum 
winds 24 hours in advance using only data which arc measured by satellites, but assuming 
that the forecaster knows whether the stprm has maximum winds >65 knots. This can 
quite reliably be determined by other techniques which use only satellite imagery (Dvorak, 
1973). The mean Tgg was used for reasons already advanced and its variability also con- 
tributed. The degree to which the lower. te^eratures are evenly distributed around the 
storm center in the mean illustrate organization and should be a measure of the efficiency 



of Uie heat engine of the hurricane to convert lieat energy into kinetic energy, This 
factor was taken into account somewhat by using averaged around a ring rather than 
for individual sectors. Another measure of distribution is the standard deviation of the 
Tgg in the ring. Even for intense hurricanes there will be areas of weak or no convection 
between cloud bands (hence relatively high temperatures). The important point is that wo 
measure how well the cold temperatures (i.e., the active convective cells) are distributed 
throughout the various octants. Therefore, we used the standard deviation of the mean 
sector temperatures os another parameter to distinguish between intensifying and weakening 
tropical cyclones. 

For the weaker storms (Vq < 65kts), the equation developed was 

AV+24w = 167.16 - 0.682T2,3 (3) 

v/here V+ 24 ^is the predicted change in maximum winds in the next 24 hours. 

For the intense storms (Vq > 65 knots), the following two equations were tested. 

V,24, = 378.51 - 1.225Ti,2,3 W) 

V+24, = 390.72 - 1.246Ti^ 2,3 - 0.506a3 (5) 

where 03 is the standard deviation of the T^g in the eight sectors of ring 3. The mean 
errors for forecasts by these equations for both dependent and independent data are given 
in Table 2,'*' The verifications listed in Table 2 for Equation (3) show it should be quite 
useful for identifying the weak storms that are intensifying. Results are significantly better 


*Table 2 includes verifications of several equations wlUch vrerc tested including the information contained 
in Table 1 . The latter was repeated to facilitate comparison of results. 
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Tabic 2 

Mean Errors of Forecasts 


Dependent Data 
(1970 Storms) 

Indep 
(1973 an 

endent Da 
d 1974 St 

tu 

orms) 

Erpiation 

n 

E 

E-NC 

E-P 

11 

E 

E-NC 

E-P 




Weak Tropical Cyclones 







(Vo < 65 Knots) 




1 

Hj| 

8.6 

-4.0 

-4.2 

|M| 

12.6 

-4.2* 


3 

KM 

9.6 

-3.0 

-3.2 


12.8 

-4.0* 


6 

24 

24.9 

12.3 

12.1 

33 

25.0 

+8,2 

+11.8 

7 

24 

10,2 

-2.4 

-2.6 

33 

14.2 

-2.4 

mm 




Intense Tropical Cyclones 







(Vo > 65 Knots) 




2 

34 

9.4 

-7.5 

-2.6 

o 

16.1 

-3 1 

-6.9* 

2(-llkts) 





Bl 

11.8 

-7.4* 

-11.2* 

4 

34 

11.4 

-5.5 

-0.6 

20 

18.8 

2.0 

5.6 

5 

34 

10.8 

-6.1 

-1.2 

20 

18.7 

1.9 

5.5 

6 

34 

21.1 

4.2 

9.1 

20 

17.4 

-1.8 

-5.6 

7 

34 

12.1 

-4.8 

0.1 

20 

14.7 

-4.5* 

-8.3** 

4 (-1 1 kts) 





20 

15.6 

-3.6 

-7.4* 

5(-llkts) 





20 

14.9 

-4.3 

-8.1* 

All Storms 

1 &2 

58 

IQI 

-6.0 


53 

13.9 

-3.8* 

-3.0** 

6 

58 


7.6 


53 

22.1 

4.4 

+5.2 

7 

58 

11.4 

-3.7 

-0.9 

53 

14.4 

-3.3* 

-2.5 


NOTES: E is mean error of forecasts made by indicated regression equation. 


E-NC is mean difference in errors of forecasts made by regression equation and 
by assuming no change in wind speed during forecast period, 

E-P is mean difference in errors of forecasts made by regression equation and 
by assuming the cliange in wind speed in the next 24 hours would be the same 
as the change during the previous 24 hours. 

♦Differences significant at 1% level. 

** Differences significant at 5% level 
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at the 1 percent level than forecasts of no change and superior to forecasts of persistence 
of change for the independent scries. Equation (4) made good forecasts for tlie t. -pendent 
data but tlicy were inferior to those made by the persistence forecasts for the independent 
data series, Tliis was also true for Equation (S), but a variation of this equation which wili 
be discussed in Section 6 made forecasts much superior to tliosc made by tlic persistence 
techniques even in the tests witli independent data. The great value of these forecast 
equations is ttiat no data ottier tiian the temperatures observed by satellite are needed ex- 
cept for the approximate location of the center. The center can be selected by inspection 
of the imagery and the grid may then be placed to compute mean temperature values. 

Since no knowledge is required of the previous or current intensity of the storm, except 
for intensity classification (over/under 65 knots), verification given in Table 2 indicates 
that Equatifms (3) and (5) could be very useful in areas where only satellite data are 
available. 

Equations were next developed and tested which apply to all the storms, that is, they 
do not require classification of the storms into weak and intense categories. The accuracies 
(see Table 2) are good enough to suggest the equations can be very useful, especially in 
areas where only satellite data are readily available. These equations are: 

V+24^ = 200.51 - 2.213ai_2 - 0.381 Tj , 3 ( 6 ) 

V+24a = 146.6 - 1.569aj^2 + 855Vo - 0.513Ti (7) 

In tests on the independent data the mean of the errors is 22.1 knots for Equation (6) 
which uses only temperatures and their distribution, but is only 1 4.4 knots for Equation 
(7). Tlie latter is significantly better at the 1 percent level than those of forecasts of no 
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cliiingc ami (s superior to thosi? of forccnMs of persistence of change (13 percent level of 
significance). 

Equation (7) makes tlie best forecasts for tlie intense storms of tlie indepciulcnt scries 
of any of the eqimt^ .presented in this paper, nitliougli forecasts for Equation (2) were 
slightly better for v .‘,0 dependent data. For tiie dependent scries of intense storms mean 
error for Equation (7) was 12,1 knots which compares with errors of 16,9 and 12,0 knots 
respectively for forecasts of no cliange and persistence of change. 

6.0 BIAS IN RESULTS 

Tlie results from Equation (2) witliout any adjustment are significantly better than 
those for persistence. However, examination of tlie errors reveals that the regression equa- 
tion forecasts maximuni wind speeds 11.3 knots too high in the mean for the independent 
cases from tlie 1973 and 1974 seasons, In this and latter sections a value, which will be 
referred to as the bias, Is calculated by taking ■':.<!! algebraic mean of the errors, If tlie 
constant in Equation (2) were reduced by 1 1 knots tlie bias of the forecasts for independ- 
ent data would be about zero and the mean error of the forecasts us defined in footnote I 
of Table 1 would be 1 1.8 knots which Is much superior to the forecasts made by assuming 
either no clumge or persistence of cliange. In fact, it is probably better tlian tlie fore- 
casters arc currently doing,* Tlie magnitude of the mean error is not especially sensitive 
to the particular value used in reducing the constant term. For example, the mean error 
varies from 12,7 to 1 1,6 knots if the constant is reduced by various values ranging from 

•It Is considerably smaller than 13,5 knots wlilch is the average error tlie forecasters had for the 1973 and 
1974 seasons. It has not been practical, however, to make comparison with what tlie forecasters did in 
tiic particular 20 cases used in tliis study. 
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8 lo 17 knots, The critical question Is then, wity iloes the equation forecast consistently 
too high values? Did the sensor calibration cimnge or was there really tiiat much change 
Iti tlie relationship between mean temperatures anil maximum wind speeds from 1970 to 
the 1973-1974 seasons? 

Since two satellites were used in collecting the dependent and independent data, 
checks were made to determine If calibration ol* the sensors varied, The Tjjj, near botli 
the cold and warm limits of tlie scale were compared to determine If there was variability 
in the observed values between the sensoni of NIMBUS IV and V. The mean Tjjy of the 
coldest sectors of the typhoon served as values for lower limit teiuperaturcs and the clear 
area sea surface temperature values functioned ns the upper limit. Tlie differences between 
the values of the two satellites were slight (<2°C) in both magnitude and ninge. Tiie data 
from one satellite were therefore considered to be cornpatible to those of tlie otlier. 

In earlier paragraphs it was explained that the relationship between storm intensity 
and mean cloud top temperatures probably varied witli latitude, season, sea surface tem* 
peraturcs and ludght of the tropopause, as well as the parameters evaluated. Table 3 com- 

M 

pares the distribution of the dependent and independent cases both temporally and spa- 
tially. It is obvious that tliere arc large differences. For example in the dependent cases 
for intense storms 88 percent of tlie cases occur in the months of August througli October, 
By contrast for the independent data only 45 percent of the cases come during tliose 
montlis. For tlie weak storms only 12.5 percent of the cases come during May, June and 
July, but for the independent data SI. 5 percent of the cases occur during those months. 
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Ijiblc 3 

Distribution of Storms by Montlis 


Montli 

Weak Storms (Vo < 65 Kt) 


Inlonso Storms {V„- 

^o5Kt) 

Dependent Data 

Independent Data 

Dependent Data 

Independent Data 

n 

Mean Error 
of Forecasts 

n 

Mean Error 
of Forecasts 

It 

Meaiv Error 
of Forecasts 

n 

Mean Error 
of Forecasts 

May, June, July 

3 

5.2 


0.7 

2 

1.2 

7 

10.7 

August 

6 

0.2 

2 

-7,8 

14 

-2.0 

3 

5.9 

September 

5 

3.2 

7 

-14,6 

7 

2.4 

3 

5.7 

October 

3 

-20.4 

6 

0.2 

9 

3.3 

3 

13.0 

November 

7 

+4.0 

1 

10.7 

2 

-10.7 

4 

19.5 

All Cases 

24 

0.0 

33 

-3.0 

34 

0.0 

20 

11.3 


Distribution of Storms by Latitude 


Latitude 

Weak Storms (Vo? 65 Kt) 

Intense Storms (V^ > 65 Kt) 

Dependent Data 

Independent Data 

Do 

oendent Data 

independent Data 

n 

Mean Error 
of Forecasts 

n 

Mean Error 
of Foiccasts 

n 

Mean Error 
of Forecasts 

n 

Mean Error 
of Forecasts 

<10® 





2 

■H 



10 < LL?20® 

19 

-0.9 

24 

-0.9 

11 


8 

11.8 

20® <LL<30® 

4 


6 

-6.6 

12 


10 

10.1 

>30® 

1 


3 

-11.8 


-D.’J 

2 

15.8 

Total 

24 


33 

-3.0 


0.0 

20 

11.3 


Distribution of Storms by Arens 


1 

Weak Storms (V^, < 65 Kt) 

Intense Storms (Vq > 65 Kt) 


Dependent Data 

Independent Data 

Dependent Data 

Independent Data 


n 

Mean Error 
of Forecasts 

n 

Mean Error 
of Forecasts 

n 

Mean Error 
of Forecasts 

n 

Mean Error 
of Forecasts 

Latitude > 1 5“ 
LonBltude>125® 

7 

-1.0 

D 

-12.4 

20 

-0.4 

14 

3.0 

Longitude <121° 


“0.4 

14 

0.8 

6 

5.9 

3 

19.3 

Latitude <28® 
Longitude > 123° 

'■ 

-0.5 

15 

0.1 

17 

-1.6 

12 

12.6 


Mean Errors (Mean Bias) 


Weak Storms (Vq? 65 Kt) 

Intense Storms (Vq > 65 Kt) 

Dependent 

independent 

Dependent 

Inde 

pendent 

IQIII 

No, 

Cli 

mi|]i 

li 

Reg. 

No 

Ch 

Pers, 

ofCli 

n 

Reg. 

No 

Cli 

Pers, 
of Cli 

n 

Reg. 

No 

Cli 

Pers. 

ofCh 

0.0 

5.1 

1.5 

33 

Bl 

-13.1 

0.5 

34 

0.0 

-1.4 

10.6 

Q| 

11.3 

0.8 

16.7 


Reg, is tiie appropriulc regression Equullon. 
No. Gi is no change, 

Pars, of Gi is Persistance of Change. 
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Tlio distribution by iutitiulc Is more comparable, iievertiieiess tberc alVtltlTcronccs, 

Tor cxaniplc witli tiie intense storms, 2 from tlie dependent series occurred at latitudes 
considerably less than 10 degrees and 9 were at latitudes greater than 30 degrees. From 
the Independent sample none were south of 10 degrees and only 2 wore north of 30 
degrees. Tliese differences between the siuhplcs motivated invcstijvUion of wlietiier lati- 
tudinal or seasonal differences in the relationship between mean temperature and maximum 
winds coidd account for the bias in the forecast errors for the 1973 and 1974 seasons, 


Tlicre are certainly differences. For the weak storms in November the bias of the 
forecasts is +4 knots for the dependent scries and +10.7 knots for the independent series. 
Likewise in October for the intense storms the bias for the dependent series is +3.3 knots 
and 13.0 knots for the independent series. For the latitude belt between 20 and 30 
degrees north latitude for the intense sloritis the bias for the dependent cases is 4.3 knots 
aitd 10.1 knots for the independent cases. Tlieso suggest that the bias in the distribution 
; might account for part of the bias in the forecast errors. There arc, however, ii lot of in- 
consistencies between the dependent anj -Independent data. For example the ni'ean error 
for the intense storms is -10.7 knots in November for the dependent and +19.5 knots for 
the independent cases. Likewise +2.3 knots in the latitude bolt of 20-30 degrees for the 
dependent cases of the Weak storins and -6.6 knots for the independent scries. Even the 
case- of 4.3 knots versus 10. 1 knots quoted earlier for the belt of 20-30 degrees north 
latitude of the intense storms is inconsistent. That is, for the dependent scries the bias is 
4.3 knots greater than the bias for all the storms in that series, while the +10.1 knots 
for the independent series is 1,2 knots less than the bias for all the storms. Because of 
these inconsistencies, further analyses Were made. Regression equations were developed 
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' for both the dependent and independent series relating tlie errors of the forecasts made 
by Equation (1) or (2) to the errors made by forecasts of no change and forecasts of 
persistence of clianges. These regression equations accounted for and 69 percent of 
the variance in the errors of tiic forecasts made by Equations (I) and (2) for the dependent 
and independent series respectively, Tliat is, the errors in forecasts Jiiadc by Equations (I) 
and (2) are liighly correlated with tlie rate at which storms change intensity or acceleration 
in the rate of change of intensity during forecast periotis. Using these relationships, further 
analyses were made by month and area of the bias in errors, In sonic cases where the mean 
error had been positive, this analysis suggested the bias of the error should be negative or 
vice versa. 

In summary, the principal conclusion is that with the number of cases available it is 
very difficult to assign any large portion of the error to latitudinal or seasonal variations. 

In most cases there were enough inconsistencies between indications from the dependent 
and independent series as to make conclusions doubtful. Tlie only cases in whicli there 
is sufficient consistency arc: • 

Weak storms: 1. May, June and July forecasts are too iiigli by a small amuUiit. (-1 ) 

2. November forecasts are too high. (-2) 

3. Area >15 degrees north latitude and > 125 degrees east longitude 
forecast values are too low. (+4) 

Intense storms: I. October forecasts arc too high. (~2) 

2. South China Sea storm forecasts are loo liigli. (-2) 

Numbers in parentheses by each case arc best estimates of the correction that should be 
applied to forecasts made by Equations (I) and (2). It should be emphasized that these 
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corrections tire not supported by iulcquntc djlta and might be quite diftcreut for another 
sample. Tl « )ould be used, if at all, witjji great' caution. It is obvious, however, from 
the magnitude of the indicated correction^/and the Iqck of indications for many of the 
cases, that one cannot account for the larfe bias in the error of the hyJependcnt cases of 
the intense storms by seasonal or locational blas'in the data sample. 

Another explanation appears more likely for the bias. Atkinson and Holliday (1977) 
discussing techniques used at Guam’s Joint Typhoon Warning Center (JTWC) to convert 
obseived central pressures to estimated maximum winds for tropical Gj/clones wrote: 

“There was considerable Uncertainty involved in the existing equations and a general belief 
among JTWC forecasters that they overestimated the maximum winds. Therefore, in 1973 
a new pressure-wind relationship developed by Fujita, et al. , (1971) was adopted for oper- 
ational use. While the Fujita relationship appeared to give more realistic wind values, a 
large-scale data collection effort (described in the next section) was initiated to obtain 
sufficient information to verify or refine the existing relationships.” 

Ciptain Holliday was a forecaster In JTWC 1970-1975. Col. Atkinson was Director 
of the JTWC in 1973 and 1974. In a personal communication Captain Holliday reports 
there was a definite change in philosophy among the JTWC forecasters by early 1974 
when they were preparing the “best track” information for the 1973 season. He says 
by then the information reported by Atkinson and Holliday (1977) was largely available 
and was used in converting 700 mb heights and sea level pressures measured by reconnais- 
sjuice aircraft to maximum winds at a lower value for the “best track” presentation for 
the 1973 and 1974 seasons than they did fo' the 1970 season. He could not estimate the 


25 


mattniUulc efioctcd by the chaiiKc bul IhouithI that applying the lakahakhi (l‘>52) and 
Atkinson and liollklay ( l‘>77) cinialions to representative data might provide a realistic 
estimate. 

lakahashi’s later work suggested that the constant in his ei)uation should be changed 
when the storms were at highci latitudes. The three equations Ivlow are trom Atkinson 
and Holliday’s paper (l‘>77). i quation (8) was developed by them, Equation (d) is Irom 
lakalushi's paper (1^52), and I'quation (10) is Takahashi’s equation for higher latitudes, 
lable 4 compares the maximum winds calculated from a range ol minimum central 
pressures by the three equations. 


Table 4 

Relations Between Minimum Central l*ressure and 
Maximum Winds in Tropical ('ycloncs 


‘V 

V.n, 


v.„, 

1000 

.10 

42 

36 

ddO 

46 

(>0 

51 

480 


73 

63 

470 

72 

85 

73 

460 

83 

45 

81 

450 

44 

104 

84 

440 

103 

1 12 

46 

4.10 

113 

120 

103 

420 

122 

127 

104 

410 

130 

134 

115 

400 

138 

141 

121 

840 

146 

1 ’7 



26 










( 8 ) 


V„,| - h.7(l01() - 
V,„j » 13.4(1010 - 

Vp,j » 11. 5(1010 - P,.)"-^ (10) 

P^. i> tlu* central pressure ot the storm in millibars, and V,„ is the calculated maximum 
wind speed in the storm at the time P^ was measured. 

In the independent sample of 20 intense storms from the 1473-1^74 seasons, 2 of 
the storms were at latitudes greater than .^0 degrees. Assuming Takahashi’s Equation (10) 
for higher latitudes applies to these 2 storms and i-quation (^1 applies to the other 18 
storms, then the differences between results from I quation (8) and the combinatic n of 
Pquations (d) and (10) in calculating maximum wind speeds for intense storms are 
as follows: 


Fable 5 

Variation of Maximum Winds 
Related to Central l*ressure 



Difference in Maximum 
Winds (Knots) 

*)H0 

12 

470 

12 

460 

1 1 

450 

4 

440 

7 

430 

5 
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The mean of the maximum winds for the inclcpemleiU data cases was 88 knots. Prom 
Equation (8) tlds means a central pressure of 95S knots. Table 5 suggests a difference 
between the maxinnun v/inds for 1970 ant! those of 1973, 1974 of 10 knots. At this 
late date it is not practical to reconstruct all of the work the forecasters did in preparing 
tlic “best track" data for 1970, 1973 and 1974. The above results however suggest that 
cliangcs in procedures at the JTWC may account for much of the bias observed. 

If the change in procedures can account for the bias observed in the forecasts made 
by Equation (2) a further clieck would be provided by comparing its constant term with 
the one in a regression equation developed using the 1973, 1974 sample as dependent 
data. An equation of llic same form as Equation (2) developed using 1973 and 1974 
data as dependent data is; 

V+24, = 243.91 - + 0.385 AV.24 + 0,424 Vq (M) 

At first inspection the constant here appeals I'd considerably larger, rather than smaller, 

'''''' -i' 

than the one in Equation (2), but the other terms conceal part of the constant. For ex- 
ample, Equations (2) and (1 1) can be expressed respectively, 

V+24s = 227.86-0.761210 + {Ti ^2,3 -210)1 io.499AV„24 + 0.398t65 + (Vq - 65)1 (12) 
V+24, = 2 43.91 -0.885 [210 t(Ti ,2,3 -210)1 +385AV_24 +0.42^1[65 +(Vp _ 65)1 (13) 
If the equations are simplified and if (13) is subtracted from (12), 

Vf24s" ^24, = 8-5 + 0.125(T\4;3 - 210) + 0.1 14AV.24 - 0.026(Vq - 65) 

If this is solved for representative values, e.g., T|,2,3 = 230“K and Vp = 90 knots, then 

= V+24g - V+24, = 10,9 knots. 

The bias for Equation (2) when tested on independent data was 1 1.3 knots. 
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Forecasts muclc by Hquation (2) with the constant reduced by 11.3 and those inacic 
witli Equation (1 1) have a correlation of 0.9899 and tlic bias in tlie forecasts is -0.015. 
Tlie mean difference In the absolute value of the errors is 1 .6 knots and the standard 
deviation of these absolute values is 2.1 knots. 

Results of the testing are quite convincing that both the 1970 and 1973-74 series 

« 

show a real relationship between the Tdd around a typhoon and the change in maximum 
winds during the next 24 hours. Furthermore, the relationships for the 2 periods agree 
quite closely except for the constant term. Further evidence that the change in pro- 
cedures at J'rWC account for much of this bias is as follows; 

Checks were next made with Equation (4) which did not use the change in maximum 
winds during the preceding 24 hours as a predictor. Another equation, (14), was derived 
in the same format using 1973 and 1974 data as dependent. 


V,24s = 378.5! 

- 1.225'f,,2,3 

(4) 

V+24, = 398.9 

- 1-375Ti.2,3 

(14) 


If (4) and (14) arc treated as in the preceding paragraphs and subtracted, 

~ V+24 = ll.ll + 0.15 (Ti,2,3 ~ 210) (16) 

When T, 2,3 = 230°IC is substituted, the result is 14.1 knots. Tlie bias of forecasts made for 
1973-74 data with Equation (4) was 13.6 knots. 

Likewise, the more complete Equation (5) was tested by comparing with a similar 
equation, (17), developed using 1973-1974 data as dependent. 
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V+24, » 395.3 - 1.333T,,2,3 - 0.4,04a3 (17) 

When the equations are subtracted, 

y+24s - V+24, = ‘3.078 + 0.087(T,,2,3 - 210) - 102(03 - 6) 

whicli for Tj 2,j ~ 230®K and 03 = 10, becomes 14.4 knots, The bias of forecasts made 
for 1973-74 data with Equation (5) was 13.7 knots. 

Mean errors for forecasts inatlo with Equations (4) and (5) for the 1973-74 series 
when constants in equations were reduced by 1 1 knots were 15.6 antJ 14.9 knots re- 
spectively. Reference to Table 2 shows that both of these mean errors are considered 
smaller than those of forecasts made by assuming eitiicr no change or persistence of 
change ami they arc significantly better at the one percent level than tlic errors made by 
forecasts assuming persistence of cliange. 

Sinniiated operational tests were also made of ttie regiessfnn equations. The wind 
speed data used thus far in developing and checking the regression equations were all 
fyken from the “best track” information prepared at JTWC. Wind tlata for tlie independ- 
ent series have been published in the Annual Reports prepared at the JTWC (U.S. Fleet 
Weatlier Central/JTWC, 1970, 1973 and 1974) not only in the “best tiack” reports but 
also as observed on reconnaissance fliglits and as used in the typhoon advisories prepared 
at the time. Times of these latter reports did not coincide exactly with time of satellite 
temperature data. The advisory prepared within a few hours (less than 6 hours in all 
ciises) after a satellite pass contained an estimate of the maximum wind based on informa- 
tion collected at approximately the same time as the satellite overflight. The data for 

V( 
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Vq, AV. 2 .J tuul AV^ 24 > therefore, were taken from the advisories and used to test Equations 
(I ) and (7) for the weak storms and Equations (2) and (7) for tlio storms with Vq > 65 
knots. In all eases the differences between the mean errors of the regression equations and 
corresponding mean errors of persistence forcaists were greater than those listed in Table 2. 
That is using simulated operational data caused greater increases in the errors of persistence 
forecasts than for the regressions equations. 

The forecasts made at JWC for the sjime days as the forecasts made by the regression 
equations for the independent data series are not readily compared with the latter. Their 
starting times coincide with tlie advisory times which usually differ from the times of the 
satellite passes. Furthermore, the corresponding official forecast was not always included 
in the annual reports of JT\VC. To the extent information was available, however, wc 
verified that the mean error of the official JTWC forecasts for the storms having Vq > 65 
knots for the days used in the experimental forecasts reported in this paper were much 
larger than the mean errors of tlie JTWC forecasts for the entire 1973 and 1974 seasons. 
This helps explain why the mean errors of alt the forecasts (by regression equation, per- 
sistence and no change) were larger for the independent data than they were for the de- 
pendent series. That is, the independent scries was biased toward situations which were 
difficult to forecast by procedures normally used. 

7.0 CONCLUSIONS 

Evidence presented in preceding paragraphs strongly suggest that the relationship 
between satellite measured T^jj around tropical cyclones and changes in maximum winds 
of the storms during the next 24 hours is real and is sufficiently strong to be a useful 
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forecast tool. TgQ can be used to obtain reasonably good results when used alone, and 
can greatly enhance the value of other indications of intensity changes when other signif- 
icimt data are available for use by the forecasters. 

More work, however, needs to be done to evaluate at least the constant in the equa- 
tions. Besides the variation between years for the Western North Pacific data which is 
probably accounted for at least in part by a change in procedures of the forecasters, there 
may well be a variation between oceans. The data presented in Table 3 suggest that there 
is at least a small variation by months and by geographical location. It certainly would 
seem logical that there would be some variation, for the Tub of the cloud tops are probably 
not only a function of the strength of tlie convection, but also a function of the sea tem- 
perature, the lapse rate in the troposphere, and the height of the tropopause. All of these 
Vary in the mean With the season, and with the geographical area. Certainly the proximity 
of land could also have some influence on the relationship. 

The results presented herein are sufficiently promising that it would seem desirable 
to further investigate the relationship between satellite data and storm intensity with more 
cases which would permit further stratification of the data. In particular, data need to be 
examined from the principal tropical cyclone regions of the other oceans in the world to 
see if there is any great variability between regions. Some data collected from the Western 

Atlantic and Eastern Pacific suggest that the is a useful forecast parameter, but too 

>' ■■ 

few cases from those regions have been examined to state stronger conclusions. 
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FIGURE CAPTIONS 


Figure 1, Grid used in amilyzlng the temperature data. Tlie concentric rings are spaced 
1 1 1 km apart. Tljc center of this grid coincides witli the center of the tropical 
cyclone, 

Figure 2. Comparison of mean temperatures of cloud tops around 16 imrrlcancs \vltl» 
mean temperatures from 19 storms of less tlian hurricane intensity. The 
tropical cyclones all occurred in the Western North Atlantic in 1969, The AT 
values are positive when the temperatures in the hurricanes arc liighcr, 

Figure 3. Comparison of mean temperatures of cloud tops around 1970 tropical cyclones 
of the Western North Pacific: 1 5 storms with maximum winds less than 65 
knots (broken line); 13 typhoons with maximum winds less than 100 knots 
(tong“short-long dashes); and 14 typhoons with maximum winds equal or 
greater than 100 knots (solid Mne). The insert contains graphs of the differ- 
ences in temperatures of the latter and the two weaker categories. Tlie rings 
arc illustrated in Figure 1. 

Figure 4. Same as Figure 3 except that the data are stratified according to change of 

4 

# 

intensity during next 24 hours: Maximum winds increasing 10 or more knots 
(I), maximum winds changing less than 10 knots (S), maximum winds de- 
creasing 10 or more knots (W), and storms which never reached hurricane 
intensity (T.S.). Only storms located south of 30°N were included. The in- 
sert shows that the intensifying storms have much colder cloud tops within 
8 degrees (888 km) of the storm center than the others. 


FIGURE CAPTIONS (Cont!wucd) 


Figure 5. Temporal changes of mean equivalent blackbody temperatures (rings 2, 3 and 
4 from Figure 1) and maximum winds of Typhoons Hope (Sept. 1970) and 
Billie (August 1970). Tile maximum wind changes lag the temperature changes 
in both cases, Tlie temperature scale is inverted to show the lower tempera- 
tures at the top of the graph. 

Figure (5. Comparison of Vertical Motion and Central Pressure in a Model Hurricane 
Sliowing Lag with Time (Rosenthal 1978) 



TABLE CAPTIONS 


Tabic 1 Mean Errors of Forecasts (Knots) 

Table 2 Mean Errors of Forecasts 
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